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We consider the critical point of two mean-held disordered models : (i) the Random Energy 
Model (REM), introduced by Derrida as a mean-field spin-glass model of N spins (ii) the Directed 
Polymer of length N on a Cayley Tree (DPCT) with random bond energies. Both models are known 
to exhibit a freezing transition between a high temperature phase where the entropy is extensive 
and a low-temperature phase of finite entropy, where the weight statistics coincides with the weight 
statistics of Levy sums with index n = T/T c < 1. In this paper, we study the weight statistics at 
criticality via the entropy S = — Wj In Wj and the generalized moments Yu = ^2 w i > where the Wi 
are the Boltzmann weights of the 2^ configurations. In the REM, we find that the critical weight 
statistics is governed by the finite-size exponent v = 2 : the entropy scales as Sn(T c ) ~ N 1 ^ 2 , 

the typical values e lnYk decay as N~ k ^ 2 , and the disorder-averaged values Yk are governed by rare 
events and decay as iV~ 1//2 for any k > 1. For the DPCT, we find that the entropy scales similarly 
as Sn(T c ) ~ N 1 / 2 , whereas another exponent v' = 1 governs the Yk statistics : the typical values 
e lnY k d eca y as J\f~ k , the disorder- averaged values Yk decay as iV -1 for any k > 1. As a consequence, 
the asymptotic probability distribution ¥jv =00 (g) of the overlap q, beside the delta function S(q) 
which bears the whole normalization, contains an isolated point at q = 1, as a memory of the delta 
peak (1 — T/T c )S(q — 1) of the low-temperature phase T < T c . The associated value 7Tjv=oc(q = 1) 
is finite for the DPCT, and diverges as Tf N=00 (q = 1) ~ N 1/2 for the REM. 

PACS numbers: 75.10.Nr, 02.50.-r, 05.40.Fb 



I. INTRODUCTION 

Spin-glasses 0, Q and directed polymers in random media [H are two kinds of disordered models where the relations 
between finite dimensional models and mean-field models have remained controversial over the years. For the directed 
polymer in a random medium of dimension 1 + d, with d = 3 where a localization/delocalization transition occurs, 
we have recently found numerically that the weight statistics is multifractal at criticality in close analogy with 
models of the quantum Anderson localization transition, where the multifractality of critical wavefunctions is well 
established In this paper, our aim is to study the critical weight statistics in the mean-field version of the model, 

namely the Directed Polymer on a Cayley Tree (DPCT) [7, 8]. This model presents many similarities 7, 8] with the 
Random Energy Model (REM), introduced by Derrida as a mean-field spin-glass model [9| : both models are known 
to exhibit a freezing transition between a high temperature phase where the thermodynamic observables coincide with 
their extensive annealed values, and a low-temperature phase of finite entro py, in which the weight statistics [Tol. ITlT] 
coincides with the weight statistics of Levy sums with index /i = T/T c < 1 [12j |. Therefore we also study the weight 
statistics of the REM at criticality, as well as in Levy sums at the critical value fi c = 1, to compare with the results 
for the Directed Polymer on a Cayley Tree. We find that the three models have different critical finite-size properties. 

The paper is organized as follows. In Section|TTl we recall the main properties of the Random Energy Model (REM) 
and of the Directed Polymer on a Cayley Tree (DPCT). We then study in parallel the weight statistics at criticality 
for the two models : we describe the properties of the entropy in Section [ill! the decay of disorder- averaged values 
in Section ITVl and the decay of typical values Y k VP = e lnYfc in Section [V] To explain the differences between averaged 
and typical values, we discuss in Section PVTl the probability distributions of the maximal weight w max and of Yi over 
the samples, with a special emphasis on the rare events that govern averages values. The Section IVIII is devoted to 
the finite-size properties of the overlap distribution at criticality. Section [Villi contains our conclusions. For clarity, 
the weight statistics of Levy sums is discussed separately in the Appendices : we first recall in Appendix [A] the results 
for jj, < 1, and we describe the critical case fi c = 1 in Appendix IBl 



II. MODELS AND OBSERVABLES 
A. Reminder on the Random Energy Model 

The Random Energy Model (REM), introduced in the context of spin glasses is defined by the partition function 

2 N 



Z N =Y J ^ Ei (1) 



i=l 



where the energies Ei of the 2 N configurations of N spins are assumed to be independent random variables drawn 
from the Gaussian distribution 



P N {E) = ~L=e~^ ( 2 ) 



This model presents a freezing transition at Q 



Tc = (3) 



and we now briefly recall the main properties of the high temperature and low temperature phases. 

1. High-temperature phase 

In the high temperature phase T > T c , the free-energy per spin coincides with the annealed free-energy 0] 

f(T>T c ) = f ann (T)=-Tln2-± (4) 
As a consequence, the entropy per spin vanishes linearly at the transition 

s(T > T c ) = s ann (T) = In 2 - -L « (T-T c ) (5) 
and the specific heat per spin remains finite 

1 

2T 2 T->T+ 



c(T > T c ) = c a „„(T) = ^5 _ ^_ , c(T+) = 2 In 2 (6) 



2. Low-temperature phase 

From the thermodynamic point of view, the low-temperature phase T < T c is completely frozen, with a constant 
free-energy per spin § 

f(T < T c ) = -Vha2 (7) 
As a consequence, the entropy per spin vanishes in the whole low-temperature phase 

s(T <T C ) = (8) 

as well as the specific heat per spin 

c(T < T c ) = (9) 

To understand better the properties of this low-temperature phase, it is convenient to study the statistical properties 
of the configurations weights in the partition function (Eq. [1]) 
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It turns out that their statistics is in direct correspondence with Levy sums of index < fx — ^- < 1 [l2j (see 
Appendix [A"l for more details). In particular, the moments 

2 N 
i=l 

have for disorder-averages [12] 



(11) 



r(fc)r(i- M (r)) /v ; r c 



The density /(w) giving rise to these moments 



dww k f(w) 



reads 12 



w 1 M (l — u>' 



r( M )r(i - m) 



(12) 



(13) 



(14) 



and represents the averaged number of terms of weight w. This density is non-integrable as w — > 0, because in the 
limit N — > oo, the number of terms of vanishing weights diverges. The normalization corresponds to 



F fc= i = / dwwf(w) = 1 



(15) 



In particular, as the transition is approached /i = T/T c —* 1 these disorder-averaged moments all vanish linearly 
for k > 1 



Y k cx (T c -T) 

t— >t<t 

The link between these weights properties and the thermodynamics is via the total entropy [l£ 



(16) 



Sn = - 2J Wi ln ( W i) = _ 



fc->i 



From Eq. [T2l the disorder-averaged value over the samples reads 



s n (t < t c ) = - [a fe F fe ]^ 1 = r'(i) - r^lff) 



In the critical region T — > T c , the entropy per spin presents the following finite-size scaling behavior 

— — _ ggr) 1 

SJv(T) = ^v-t"t c -^(7W) 
Similarly, the disorder-averaged specific heat per spin follows 

— CMT) 1 

CAf(J ) = — — cx 



AT t^t- N(T C - T) 2 



(17) 



(18) 



(19) 



(20) 



These finite-size scaling behaviors are in agreement with the more detailed finite-size corrections of the free-energy 
computed in Q. 
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B. Reminder on the directed polymer on a Cayley tree 

The directed polymer on a Cayley tree with disorder has been introduced in 0] as a mean-field version of the 
directed polymer in a random medium Q . The model is defined by the partition function 

Z N =Y,e~ PE{ ® (21) 
c 

where the 2 N configurations C are the paths of N steps on a Cayley tree with coordination number K = 2. The energy 
E(C) of a path is the sum of the energies of the visited bonds. Each bond has a random energy drawn independently, 
for instance with the Gaussian distribution 

p( e ) - (22) 



As in Eq. \TU\ it is convenient to consider the configurations weights w(C) — e ^ E ^ c '/Zpi in the partition function 
(Eq. l2Tj) and the associated moments Yfc (Eq. [TTj) . 



1. Similarities with the Random Energy Model 

This model presents many similarities 0, HI with the Random Energy Model described above. It presents a freezing 
transition at 

T c = ; 1 (23) 
The free energy per step coincides with the annealed free energy above T c and is completely frozen below 

/CO = / ann (T) = -Tln2-^ for T>T C (24) 



f(T) = -V21n2 for T < T c (25) 

So at the thermodynamic level, all properties are the same as in the REM : the entropy per step vanishes linearly at 
the transition (Eq. O, and the specific heat per step presents a jump (Eq. [5]). From the finite-size behavior of the 
free-energy for T < T c (Eq. 76 of 0]), one obtains by differentiation with respect to the temperature that the entropy 
per step and the specific heat per step have the same finite-size scaling as in the REM (Eqs and [201 ) 



_ S N (T) 1 



s N (T) (26) 



C N (T) 1 



Cn{T) = — t ^-n(t c ~ty (27) 

It turns out that even beside these thermodynamic quantities, the two models are still very similar. In [3], it was 
shown that, in the limit N — > oo, the distribution of the overlap q between two walks on the same disordered tree is 
simply the sum of two delta peaks at q — and q = 1 in the whole low-temperature phase [?J : 

n(q) = (l-Y 2 )6(q) + Y 2 8(q-l) (28) 

And the distribution of Yj> over the samples is exactly the same as in the Random Energy Model [7[. In particular, 
its averaged value reads 



Y 2 (T) = l- ¥ (29) 



2. Differences with the Random Energy Model 



As explained in details in [f|, the differences with the Random Energy Model (REM) is that the directed polymer 
on a Cayley tree corresponds to a Generalized Random Energy Model (GREM) with p = N levels, whereas the REM 
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corresponds to the case of a GREM with p = 1 level. This induces some differences for the finite-size properties of 
the free-energy Q . The conclusion of [|[ is that the finite-size scaling behavior of the REM only involves the product 
(T c - T)^/" with 

u = 2 (30) 

whereas for the Cayley tree, the situation is more subtle, and the product (T c — T^iV 1 /^ with another exponent 

v' = 1 (31) 

also appears. As a consequence, even if the weight statistics is the same in the low-temperature phase of the two 
models, their critical properties might be different, as we will indeed find in the following. 

C. Numerical details 

The numerical results given in the following sections for the REM and the DPCT have been obtained for the 
following sizes N (with 2 N configurations) and the corresponding numbers n s (N) of disordered samples 

N = 5 - 14,16,18,20 (32) 
n s (N) = 10 7 ,4- 10 6 ,10 6 ,4- 10 5 (33) 

III. ENTROPY AT CRITIC ALITY 

A. Disorder averaged entropy at criticality 

We have recalled in the previous Section that both for the Random Energy Model and for the directed polymer on 
the Cayley tree, the freezing transition corresponds to a jump in the intensive specific heat 



c N (T>T c ) ~ Constant (34) 

1 



c N (T<T c ) ~ — - (35) 

t— t- N(T C - Ty 

The finite-size scaling thus involves the exponent v = 2 

cjv(T) = C ((T c - T)N 1/lJ ^ with v = 2 (36) 

For the REM, the explicit form of the scaling function can be obtained from Eq. (64) of @. 

Similarly, the total disorder-averaged entropy has the following behaviors on both sides of the transition 



S N {T>T C ) ~ N{T~T C ) (37) 

1 



S N (T<T C ) ~ — — (38) 

T^T~ J-c — J 

One thus expect the following finite-size scaling form 

S N (T) = N^ 2 S ((T c - T)^ 1 ^) with v = 2 (39) 

At criticality, we thus expect both for the REM and for the Directed Polymer on the Cayley tree 

S^TJ cx N 1 ' 2 (40) 
in agreement with our numerical simulations for both models. 
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FIG. 1: (Color on line) REM : entropy probability distribution P N (S) for sizes N = 5,6, 7, 8, 9, 10, 12, 14, 16, 18, 20 (a) at T c , 
the distribution remains broad around the averaged value S oc N 1 ^ 2 , with a slow decay of rare events of small entropy S ~ 0. 
(b) at T = 2. > T c : the width around the average value S = Ns ann (T) decays exponentially in N, in agreement with Ref. [l4| ] 




FIG. 2: (Color on line) DPCT : entropy probability distribution P N (S) for sizes N = 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20 (a) at 
T c , the distribution remains broad around the averaged value S oc JV 1//2 . The statistics of rare events in the region S ~ is 
different from the REM (see Fig. [T]a ). (b) at T = 2. > T c : the width around the average value S = Ns a nn(T) converges 
towards a constant, in contrast with the REM (see Fig. Q]b ). 



B. Entropy distribution at criticality 



We shown on Fig. Q] and [2] the probability distribution Pn(S) of the entropy S for the REM and for the DPCT 
respectively, both at criticality and in the high-temperature phase for comparison. At criticality, Pn(S) remains broad 
around the averaged value S oc TV 1 / 2 , with a slow decay of rare events of small entropy S ~ 0. The comparison of 
FigsQ] (a) and [2] (a) show that these rare finite samples that are still 'frozen' at T c do not obey the same statistics in 
the REM and in the DPCT (see the more detailed discussion on rare events in section IVI C|) . In the high temperature 
phase, the width around the average value S = Ns ann (T) decays exponentially in N in the REM [14], as shown on 
Fig. Q] (b), whereas it converges towards a constant in the DPCT, as shown on Fig. [2] (b). 
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IV. DECAY OF DISORDER AVERAGED VALUES Y k (N) AT CRITIC ALITY 

A. Special case k = 2 

As already mentioned in Eq. [29] for k = 2, the explicit expression of Y^iV) is particularly simple in the low- 
temperature phase, 



Y 2 (T<T C ) = 



T c -T 
T, 



In the REM where the only finite-size scaling exponent is v = 2 (Eq. [30]) , one thus expects at criticality 

1 



Y 2 (T C ,N) 



3C 



REM 



(41) 



(42) 



in agreement with our numerical simulations, as shown on Fig. [3] For the DPCT however, we find numerically that 
the decay of Y% is governed by the exponent v' = 1 (Eq. [31]) at criticality 



as shown on Fig. [3] 



Y 2 (T C ,N) 



1 

cx — 

DPCT N-*oo N 



(43) 



In Y 



-1.5 



-2.5 - 




DPCT 



1.5 



2.5 



FIG. 3: (Color on line) At criticality, the slope of In ^(iV, T c ) as a function of (In TV) (here 5 < N < 20) is of order 1/v = 0.5 
for the REM (□) and of order 1/v' = 1 for the DPCT (Q) 



B. Other values of k 



For arbitrary k, the explicit value (Eq. [T2]) can be expanded in (T c — T)/T c as follows 

(T c - T) 



>fc(T < T c 



[l + 0((T c -T))\ 



(k - l)T c 

For the REM where the only finite-size scaling exponent is v = 2 (Eq. [30]) . one thus expects at criticality 



REM (k - l)jVV2 



1 + 



1 



ATi/2 



(44) 



(45) 



in agreement with our numerical simulations. For the DPCT, we find numerically that it is the exponent v 1 = 1 (Eq. 
I3T]) that governs the critical behavior 



Y k {T Cl N) 



DPCT N 



l + O 



JV 



(46) 
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V. DECAY OF TYPICAL VALUES Y* yp (N) = e lnYk AT CRITIC ALITY 



From the explicit expression of Eq. IA19I of lnYfc in the low-temperature phase with \x = T/T c , one obtains the 
following expansion in (T c — T) 



lnr fc = k [1 + oi(T c —T) + 0((T C - T) 2 ] ln(T c - T) + b + h(T c - T) + 0((T C - T) 2 ) 
For the REM where the only finite-size scaling exponent is v — 2 (Eq. [50]) . one thus expects at criticality 



lnY k {T c ,N) 



k 

oc — — 

REM 2 



1 + 



In N + est + O 



in agreement with our numerical simulations. 

For the DPCT, we find that it is the exponent v' = 1 that governs the decay of typical weights 



In Y k (T c ,N) 



\nN + cst + { — 

N 



DPCT 

To summarize, the exponents governing the decay of typical values are exactly linear in k in both models 

1 



[Y^(N)] 



REM 



iV fc /2 



and 



[y k tvp {N)} 



= e lnYk cx — 
DPCT c ^ jyfc 



(47) 



(48) 



(49) 



(50) 



(51) 



in contrast with disorder averaged values, where the exponents do not depend on k fEqs [431 and 14^1). To explain this 
difference, we now discuss the histograms of w max and of Yi at criticality. 



VI. PROBABILITY DISTRIBUTIONS OF uw* AND OF Y 2 AT CRITICALITY 



A. Probability distribution of w max {N) at criticality 



For each sample, we consider the maximal weight 

Wmax = maxi [wi] (52) 

among the 2^ configurations. We show on Fig. [4] and [5] the probability distribution P/v (In w max ) over the samples 
for the REM and the DPCT, both at criticality and in the high temperature phase for comparison. At criticality, 
P/v(lnw maa: ) remains broad around the averaged value 



lnw max ~ — (lnN) + ... for the REM (53) 

In Wmax ~-(lniV) + ... for the DPCT (54) 

with a slow decay of rare events near the origin hiw max ~ 0. Again, as for the entropy distribution (see Figs[T] (a) 
and [5] (a) ), the statistics of these rare 'still frozen' samples is not the same in the REM and in the DPCT as shown 
on Figs. S](a) and[5](a). In the high temperature phase, the width around the average value ln-W max oc —N converges 
towards a constant in both models, as shown on Fig. [4] (b) and [5] (b). 
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FIG. 4: (Color on line) REM : Probability distribution Pjv(lnw mal ) of the maximal weight among the 2^ configurations, 
for sizes N = 5,6, 7, 8, 9, 10, 12, 14, 16, 18, 20 (a) at T c : the distribution remains broad around the averaged value lnui maa; ~ 
— (ln7V)/2, with a slow decay of rare events near the origin \nw ma x ~ 0. (b) at T = 2 > T c : the width around the average 
value \nw max oc ~N converges towards a constant. 




FIG. 5: (Color on line) DPCT : Probability distribution Pjv (In Wmaz ) of the maximal weight among the 2 N configurations, 
for sizes N = 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20 (a) at T c : the distribution remains broad around the averaged value \nw max ~ 
— (In TV), with a slow decay of rare events near the origin In Wmax ~ 0. (b) at T = 2 > T c : the width around the average value 
In w m ax oc — N converges towards a constant. 



B. Probability distribution of Yi at criticality 



We show on Fig. [6] and [7] the probability distribution Pjv(ln Y2) over the samples for the REM and the DPCT, both 
at criticality and in the high temperature phase for comparison. At criticality, (In Y 2 ) remains broad around the 
averaged value 

hTTa ~-(lniV) + ... for the REM (55) 
\nY2 ~ -2(lnJV) + ... for the DPCT (56) 

(57) 

with again a different decay of rare events near the origin In Y2 ~ 0. In the high temperature phase, the width around 
the average value InY^ oc — N converges towards zero for the REM, as shown on Fig. [H] (b), and towards a finite 
constant for the DPCT as shown on Fig. [7] (b). 

As a final remark, it is interesting to compare the probability distribution of lnY" 2 at criticality for the directed 
polymer on the Cayley tree (Fig. [7] a) and in dimension 1 + 3 (see Fig. 3 a of [4]), where as N grows, the distribution 
simply shifts along the x-axis with a fixed shape. 
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FIG. 6: (Color on line) REM : Probability distribution P N (ln Y 2 ) for sizes TV = 5,6,7,8,9,10,12,14,16,18,20 (a) at T c : 
the distribution remains broad around the averaged value \sxY2 — —(In TV), with a slow decay of rare events near the origin 
In Yi ~ 0. (b) at T = 2 > T c : the width around the average value I11I2 c< —N converges towards zero. 
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FIG. 7: (Color on line) DPCT : Probability distribution Pjv(m y 2 ) for sizes TV = 5,6,7,8,9,10,12,14,16,18,20 (a) at T c : 
the distribution remains broad around the averaged value lnY^ — — 2(lnTV), with a slow decay of rare events near the origin 
In Y2 ~ 0. (b) at T — 2 > T c : the width around the average value In Y2 oc TV converges towards a constant. 



Rare events where w max ~ 1 and Y% ~ 1 



In the low temperature phase T <T C , the statistical properties of w max and Y2 have been studied in details in 
In particular, the probability distribution PT<T c (w m ax) coincides for 1/2 < w max < 1 with the weight density f(w) 
given in Eq. 1141 Near the transition fj, — T/T c — > 1 _ , the singularity near w max reads 



T<T, 



{Wmax) = fT<T c (w m ax) 



(1 - - W max ) 



In this section, we are interested in the behavior of the probability distribution PT c ,N(wmax) near w r , 
samples at criticality 



PT c ,N(w ma x) = fT a ,N{w m ax) - 

The same singularity governs the probability distribution of Y 2 



Aff(l - WmaxY 



Pt c M y z) - A N {l-Y 2 ) 
r 2 -»i 



(58) 
1 for finite 

(59) 
(60) 



The amplitude An represents the global scaling of the rare samples which are 'still frozen', whereas the exponent 
a describes the shape of the singularity. These rare events govern the disorder-averaged values Yfc at criticality (Eq. 
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FIG. 8: (Color on line) REM at criticality : statistics of rare finite samples which are still 'frozen' at T c for sizes N = 
5,6,7,8,9,10,12,14,16,18,20 (a) Probability distribution of the maximal weight w m ax in the region 1/2 < w max < 1 (b) 
Probability distribution of Y2 in the region 1/2 < I2 < 1 
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FIG. 9: (Color on line) DPCT at criticality : statistics of rare finite samples which are still 'frozen' at T c for sizes N = 
5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20 (a) Probability distribution of the maximal weight for sizes w m ax in the region 1/2 < Wmax < 1 : 
note the difference with the REM (Fig. [8] a ) (b) Probability distribution of Y2 in the region 1/2 < Y2 < 1 : note the difference 
with the REM (Fig. |S]b) 



I13[) . and for large k, the exponent a governs the power-law dependence in k 

For the REM, where all finite-size scaling properties involve the factor (T c — T)^ 1 / 2 , we expect 



1 

CTREM = (63) 



A R N EM „« (62) 



This is in agreement via Eq. [61] with the leading behavior of the disorder-averaged values Yfe of Eq. [45] We show on 
Fig. [5] the behavior of the probability distributions of w ma x and I2 near w max — > 1 and Y2 — > 1 

For the DPCT, the situation is more subtle. From the behavior in TV of disorder-averaged values of Eq. 0^1 we 
conclude that the amplitude is governed by v' = 1 



(64) 
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However in contrast with the REM, the behavior of the probability distributions of w max and Yi near w max — > 1 and 
y 2 — * 1 shown of Fig. [H] corresponds to a value a > 1 for the singularity exponent. The measure of the fc-dependence 
of Eq. [61] indeed leads to a value of order 

a DPCT — 1.5 (65) 

The fact that a finite a appears at criticality for the DPCT, in contrast with the REM where a — in continuity 
with the low-temperature phase, indicates that the tree structure plays a role at criticality, in contrast with the low- 
temperature phase where the overlap distribution is concentrated on q = and q = 1 (Eq. I28|) . In the next Section, 
we describe the finite-size properties of the overlap distribution at criticality. 



VII. OVERLAP DISTRIBUTION AT CRITICALITY 



In disorder-dominated phases, the order parameter is the 'overlap' between two thermal configurations in the same 
disordered sample. In this Section, we discuss in detail the overlap distribution for the DPCT, and compare with the 
REM case in the end. 

For the DPCT, we consider the probability P/v(i) that two walks of N steps have t common bonds in a fixed sample 
of a Cayley tree, where the possible values are t = 0, 1, ..N. The normalization reads 

N 

£iV(t) = l (66) 
t=o 

The usual overlap distribution ttn(q) concerning the fraction q = t/N of common bonds reads 

?r N (q)=NP N (t = Nq) (67) 

with the normalization 

/ dqw N (q) = l (68) 
Jo 



A. Reminder on the overlap distribution for T < T c for the DPCT 

As recalled in Eq. [28] the distribution of the overlap q between two walks on the same disordered tree is simply 
the sum of two delta peaks at q = and q = 1 in the whole low-temperature phase Q , and in particular the disorder 
average over the samples reads (Eq. 



7fjv=oo(9) = J%) + ( X - J) % - 1) (69) 

The finite-size corrections have been studied in [l5|,[l6[ : the probability Pjv(0 is finite at t = and at t = N, whereas 
for <C t <C N, the disorder averaged probability P/v(i) obeys the scaling 



^(0<t<iV) iv oC oo ^ 7i ^-j (70) 

where the function ip(q) presents the singularities q~ 3 / 2 and (1 — q)~ 3 ^ 2 near the two boundaries q — > and q — ► 1. 
For the finite-size overlap distribution, Eq. [7D] translates into the finite size correction 

n N (0<q<l) N ^ oo j^ m ^(q) (71) 

to the asymptotic result of Eq. [69] 



B. Finite-size overlap distribution at T c for the DPCT 

In the limit N -> 00, Eq. [69] becomes for T = T c 

7TT c ,Ar=oo(g) = S(q) (72) 

i.e. the whole normalization is concentrated on q = 0. Here we are interested into the finite-size corrections to this 
result. We show on Figs ITD1 (a) and (b) the probability distributions P Tc ^ N (t) and Tf Tc ^ N (q) for various sizes. We now 
discuss the intermediate region < q < 1 and the two limit values q = and q = 1. 
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FIG. 10: (Color on line) DPCT at criticality (a) Logarithm of the disorder averaged probability distribution Pjv(t) of the 
number t = 0,1, ..N of common bonds between two walks for sizes N = 10(O), 12(D), 14(0), 16(A), 18(>), 20(*) ( b ) test 
of the scaling form of Eq. [73] for the disorder averaged probability distribution 7fjv(?) of the overlap < q = t/N < 1 : 
hiip(q) = ln(jV ' 5 7Tjv(q)) as a function of q 

1. Region of intermediate overlap < q < 1 
For < t < 1, we find numerically that the disorder averaged probability Pjv(t) obeys the scaling 

(73) 



P JV (0«t«iV) jv oC oo ^ ViV 
or equivalently for the disorder averaged overlap distribution 7fjv(?) of Eq. [B7J 



7^(0 < q< l^oc^—Vfa) 



as shown on Fig. [TUlb. 



(74) 



2. Region of zero overlap q = at criticality 



For finite t and N — > oo, Pjv(t) converge to finite values as shown on Fig. [10] a, in particular 

P N (t = 0) ~ 0.23 
Pjv(t=l) ~ 0.15 



N^oo 



(75) 
(76) 



such that the normalization of these finite values corresponding to q = after rescaling, is 1 (Eq. [72]) . From the 
matching with the scaling regime of Eq. 173] one expects the following power-law decay for large t 



P 



N=oa(t).<X TTE 



(77) 



3. Probability of full overlap q = 1 at criticality 

By definition, the probability P/v (N) of a full overlap t — N coincides with the probability Y 2 = ^2 w 1 that the two 
walks end at the same point 



P Tc , N (t = N) = Y 2 (N,T c ) 
Using Eq. [43] the average over the samples yields 

P Tc , N (N) =%(N,T C ) oc -j= 

N^oo iv 



(78) 
(79) 
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For the disorder averaged overlap distribution of Eq. [671 we thus obtain that "Kt c ,n{(1 = 1) remains finite as N — > oo : 

W Tc , N (q = l) ~ 7f JV=oo (( Z = l)>0 (80) 



N- 



Beside the delta function S(q) which bears the whole normalization (Eq. I72[) . the asymptotic probability distribution 
T^N=oo(q) of the overlap q, thus contains an isolated point at q = 1 where TfN =oc (q = 1) > 0. This finite value at q = 1 
is due to rare events, since the typical value at q = 1 is of order (Eq. IBTj) 



4!M? = i) = ^% tw (^e) CX ^ 



typ/ 



N- 



N 



(81) 



We show on Fig. [TT]a the probability Pt c .n{^{<1 = 1)) over the samples of the probability density TT]y(q = 1) of full 
overlap between two configurations. From the probability of rare events with Y 2 — > 1 of Eq. [501 ° ne obtains via the 
change of variables ttn{<1 = 1) = -/VY2 the following singularity near the maximal value 7r(g = 1) — > A*" 



PT c ,Ar(7r( 9 = 1)) 



.4 



TT{q=l)->N N 1+<T 



[N-Tv(q=l)] c 



(82) 



where A 



(DPCT) 
N 



oc 1/N and gdpct ~ 1.5. 
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FIG. 11: (Color on line) Statistics over the samples of the probability density 7r(q = 1) of full overlap between two configurations 
at criticality, for sizes N — 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20 : (a) Logarithm of the probability distribution P(ir(q =1)) for the 
DPCT (b) Logarithm of the probability distribution P(ix{q = 1)) for the REM 



C. Overlap distribution at criticality in the REM 



In the REM with N spins, the spin overlap 

N 

^(1) C (2) 



can be defined from its relation with the p-spin glass model [12|, [l3| in the limit p — > 00. It is t = N if the two 
configurations are identical = and t < N if the two configurations are different 7^ C^ 2 \ 

As in Eq. [751 the probability Pn(N) of a full overlap t — N coincides with the probability Y2 = X) w f that the two 
configurations are the same. Using Eq. the average over the samples yields 

CW-^TJ^^ (84) 
For the overlap distribution of Eq. (67j we thus obtain that Wtc,n(.Q = 1) diverges as N — > 00 : 

W^(q = l) ^ N 1 / 2 (85) 
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Beside the delta function S(q) which bears the whole normalization (Eq. [72)) , the asymptotic probability distribution 
T^N=oo{q) of the overlap q, thus contains an isolated point at infinity as a memory of the delta peak (1 — T/T c )S(q — 1) 
of the low-temperature phase T < T c . Again this divergence at q = 1 is due to rare events. However here, in contrast 
with the directed polymer (Eq. I5T} . the typical value at q = 1 remains finite (Eq. 1ST))) 

ir% p N (q = 1) = NY* VP (N, T c ) ~ est (86) 

We show on Fig. [Tl]b the probability PT c ,N{^{q = 1)) over the samples of the probability density 7Tjv(<z = 1) of full 
overlap between two configurations. The singularity near the maximal value is given by Eq. [521 where A3 E oc 
1/N 1 / 2 and a rem = 0. 

VIII. CONCLUSION 

In this paper, we have studied the weight statistics at criticality for the Random Energy Model (REM) and for the 
Directed Polymer on a Cayley Tree (DPCT) with random bond energies. These two mean-field disordered models 
present a freezing transition with similar thermodynamic properties. In particular, between the high temperature 
phase of extensive entropy and the low-temperature phase of finite entropy, the entropy at criticality scales as Sn (T c ) ~ 
iV 1 / 2 in both models. However, the statistical properties of the weights which coincide in the low-temperature phase 
become different at the critical point. In the REM, all critical properties are governed by the finite-size exponent 
v = 2 : the typical values e lnYk decay as N~ k / 2 , and the disorder-averaged values are governed by rare events and 
decay as N~ 1 / 2 for any k > 1. In the DPCT, we find that the weight statistics is not governed by the exponent v = 2 
of the thermodynamics, but by another exponent v' = 1 that had been previously mentioned in [8j in connection with 

finite-size corrections to the free-energy below and at T c . In particular, the typical values e ln Yk decay as N~ k , and the 
disorder- averaged values decay as iV -1 for any k > 1. We have also presented numerical histograms at criticality 
for the entropy, the maximal weight w max and Y<z. We have emphasized the role of the rare samples that are still 
'frozen' at T c ( i.e. the rare samples having 5^0, w ma x ~ 1, Yi ~ 1) since it is the amplitude of these rare events 
that governs the disorder averaged values Yk as well as the overlap probability density 7Ft c ,aK<7 = 1) of full overlap 
q = 1. In particular, we have obtained that beside the delta function S(q) which bears the whole normalization, the 
disorder averaged asymptotic probability distribution 7fT c ,Ar=oo(<?) contains an isolated point at q = 1 as a memory of 
the delta peak (1 — T/T c )5(q — 1) of the low-temperature phase T < T c . The associated value TTN =00 (q = 1) is finite 
for the DPCT, and diverges as n N=00 (q = 1) ~ N 1 / 2 for the REM. 

Concerning the weight statistics at criticality for the directed polymer, let us finish by some comparison between 
the mean-field version on the Cayley tree considered here and the finite dimensional version that we have studied 
recently in We should first recall that in finite dimension d, the weights of the 0(N d ) possible spatial positions 
of the polymer end-point do not coincide with the configuration weights, in contrast with the Cayley tree where 
the end-points are in one-to-one correspondence with the 2 N configurations. In finite dimension, the probability 
distributions of the maximal weight w max and of Y% reach the values w max — 1 and Y% = 1 only for T < T gap , where 
Tgap < T c 17], whereas on the Cayley tree these two temperatures coincide T gap — T c . This is why on the Cayley 
tree, the disorder averaged values Yk for k > 1 all decay with a ^-independent exponent Yk oc 1/N representing 
the amplitude of rare events where w max ~ 1 , whereas in finite dimension, the disorder averaged values Yk decay as 
Yk oc l/jV^ -1 )- ^) where the exponents D(k) have a finite limit D(+oo) > 0. Also in finite dimension, the comparison 
with the exponents D(k) governing the decay of typical values Y, VP = e lnYfc oc l/iV^" 1 ) ^ show that the threshold 
k c between the region k < k c where they coincide D(k) = D(k) and the region k > k c where they differ D(k) > D(k) 
is of order k c ~ 2 [J], whereas on the Cayley tree, the exponents for averaged and typical values are always different 
as soon as k > 1. So the role of rare events is stronger on the Cayley tree. 

APPENDIX A: LEVY SUMS FOR n < 1 

In this Appendix, we recall some properties of Levy sums with \i < 1, since their weight statistics is the same as in 
the Random Energy Model in the low-temperature phase with \i(T) — T/T c . 
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1. Weight statistics in Levy sums 

The sum 



A I 



i=l 

of M positive independent variables (xi, ..xm) distributed with a probability distribution that decays algebraically 

A 



^)_r-_^+iI ( A 2) 



has very special property when < fj, < 1 since the first moment diverges < x >= +oo [13, [l8j| : the sum Em grows 
as M 1 /^, and the rescaled variable is distributed with a stable Levy distribution. Another important property is that 
the maximal variable x max (M) among the M variables (xi, ...xm) is also of order M 1 /^, i.e. the sum Y,m is actually 
dominated by the few biggest terms. To quantify this effect, it is convenient to introduce the weights 

m = (A3) 

and their moments 

M 

Yk = Y, w * ( A4 ) 



The link with the weight statistics in the Random Energy Model can be understood as follows. The lowest energy 
in the REM is distributed exponentially 

PextremaliE) ~ e 7 ^ (A5) 

This exponential form that corresponds to the tail of the Gumbel distribution for extreme- value statistics [l!| [2(| , 
immediately yields that the Boltzmann weight x — e~® E has a distribution that decays algebraically (Eq. IA2|) with 
exponent 

» = T 7 (A6) 

In the REM, the coefficient 7 in the exponential (Eq. IA5j) is 7 = l/T c . 

Let us also mention that in the mean-field SK model of spin-glasses, exactly the same expressions of Yfe (Eq. IA8[) 
also appear 0, HH, but with a different interpretation : the weights are those of the pure states. As a consequence, 
the parameter fi(T) which is a complicated function of the temperature vanishes at the transition n(T c ) = (only 
one pure state in the high temperature state) and grows at T is lowered towards fi(T = 0) of order 0.5 |22|. This 
is in contrast with the REM model where fi(T) — T/T c grows with the temperature from fi(T = 0) = (only one 
ground state) to [i(T c ) = 1 at the transition, where the number of important microscopic states is not finite anymore. 
Nevertheless, the expression (Eq. I A8|) for the weights of pure states means that the free-energy / of pure states in 
the SK model is distributed exponentially 



with a parameter 7 (T) = fi(T)/T. 



P(f)* (A7) 



2. Disorder-averaged moments Y^^ vv 



The averaged values in the limit M — > 00 arc finite for < [l < 1 and reads [lj 



y-Levy r(fc - fl) 

Yk - r(fc)r(i- M ) (A8) 
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Let us recall how one derives this result [12[ , since it will be useful for the critical case \x c = 1 considered in Appendix 
IB] It is convenient to exponentiate the denominator according to 12] 



(A9) 



in order to perform the average 



+oo 



dti*- 1 x k e~ tx {e- tx ) M - x 



r(fc) Jo 

For large M, the integral will be dominated by the region where t is small, and one may approximate [li 



J dxp(x)x k e" tx ~ At^- k T{k - p) 



and 



yielding 



—Levy _ MAT(k - fl) f + °° x MtfAi-Ti-n)) 



(A10) 



(All) 



(A12) 



(A13) 



leading to Eq. IA8l in the limit M — ► oo. More generally, correlations functions between Yk can also be computed [12], 
in particular 



~fp) U(i-m) +M r2(i- M ) 



(A14) 



3. 'Typical values' Y^ evy (typ) = e ln 



To compute the disorder average of the logarithm of Yj, , we first rewrite 



In Yfe = In x ij - k ln (j2 x ij 



and use the identity 



luZ= / - (e-* - e- tz ) = lim 



(i 



Using Eq IA121 we obtain 



and similarly 



M 



ME*? = 



+oo 



e -t _ I e -tx* 



M 



+ 00 /"+oo 

[i ./o 



l--)r'(l) + iln(MA(-r(-/x))) 



+ -m fM^(-r(-f)) 

|i / Lb \ rv K 



so that finally in the limit M — * oo (Eq. IA15I) 



fcjr(i-f) 



my fc = (i-* ) r'(i) + -in^_ 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 
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4. Critical behaviors near the transition point p, — ► /i c = 1 
As n — ► 1, Eq. IA19I gives the following leading term 



lnF fe ~ feln(l -n) (A20) 

i.e. the typical values vanish as 



Y k LeVV (tVP) = e WY * ~ (1 - n) k (A21) 
whereas the averaged moments of Eq. IA8I vanish linearly as 

^ = ^ + 0((l- M ) 2 ) (A22) 



as well as higher moments, for instance Eq. IA14I 

n LeVy ~^ + {{ l-»?) (A23) 



This shows that disorder-averaged values are governed by the rare events where the maximal weight w max is near 1 : 
the density f(w) of Eq. [T4l becomes for /i — > 1 

/ £eo «H (1 - riw-^il - wy- 1 (A24) 

APPENDIX B: WEIGHT STATISTICS FOR LEVY SUMS AT CRITIC ALITY fj, c = 1 

In this Appendix, we describe some results on the weight statistics for Levy sums at criticality p c = 1 to compare 
with the results given in the text for the Random Energy Model and for the Directed Polymer on a Cayley Tree. For 
He = 1, the sum Em of Eq. IA1I scales as MlnM, whereas the maximal value x max among the M variables scales as 
M [H[ : the decay of the Y k is thus expected to depend on the variable (InM). 

1. Decay of disorder averaged values Y^{M) 

We start from Eq. lATOl 



Yk = F7TT / dttk ^ xke ~ tX {e- tx ) M ~ l (Bl) 
r ( fc ) Jo 

For large M, the integral will be dominated by the region where t is small, and one may approximate 

x k e~ tx = J dxp(x)x k e~ tx ~ At 1 - k T(k - 1) (B2) 

and 

~ = \ dxp(x)e~ tx ~ e - Atln T (B3) 



e 



yielding 



2. Disorder averaged entropy 

From Eq. IBU the disorder-averaged entropy (Eq. [T7] )reads 

S(M) = -d k Y k \ k=1 = M j dt (e- tx ) M ^ J dxp(x) [(r'(l) - hit - \nx)xe~ tx ] (B5) 



Using Eq. IB31 one obtains at leading order 



S(M) cx \n(MA\nM) (B6) 
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3. Decay of typical values Y* yp (M) = e lnYk 



To compute the disorder average of the logarithm of , we start from Eqs IA15I and IA16I Using Eq IB31 we obtain 

'■+ 00 dt 



e - t _ e -MAUni j _ \ n (MA\nM) 



and as in Eq |A"T8l with fi = 1 



M 



dt 
T 



-t _ L-t^ 



A I 



o -/.■ir'ii) /,in(M^(-r(-i)) 



so that finally the leading term is (Eq. IA15[) 



( M \ ( M \ 

lnYfe = In ^2 x i ~ ^ m E x% ~ — fcln(ln 



M) 



i.e 



(lnM) fc 



(B7) 



(B8) 



(B9) 



(BIO) 



4. Finite-size scaling in the critical region 

The comparison of the results for the entropy, for the disorder-averaged values and for the typical values of the Yj. 
between the phase /i < 1 and the critical point /i c = 1 shows that the appropriate scaling variable is (1 — fi)hiM, 
corresponding to a 'finite-size exponent' 



^Levy 1 



(Bll) 



This is in contrast with the Random Energy Model where the number of configurations is M = 2 N , and the appropriate 
finite-size scaling behavior (Eq. [3D]) is (T c ~ T)N X I 2 = (1 - ^{IxiM) 1 / 2 with v REM = 2. 

5. Probability distributions of w max and of Y-j 

We show on Fig. rT2]the finite-size probability distributions of the maximal weight w max and Y% at the critical value 
fi c = 1, to compare with the corresponding figures given in the text for the REM (Fig. [8]) and for the DPCT (Fig. 

with the correspondence M — 2 N . As in Eq. [53] the behavior of the probability distribution P f _ lc .M(w ma x) near 
w m ax 1 for finite sums of M terms at the critical value [i c = 1 is of the form 

P^,M(w ma x) ^ A M (1 - w max ) a (B12) 



The amplitude Am of these rare events is the amplitude that governs the disorder averaged values Yk of Eq. [B4 

1 



A M oc 



m->oo InM 



(B13) 



The singularity exponent a is simply 



0\Le 



(B14) 



in continuity with the rare events in the region \i < 1. This value is the same as in the REM (Eq. I63[) but different 
from the value measured in the DPCT (Eq. [S5]) . 
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FIG. 12: (Color on line) Weight statistics in Levy sums of M terms ( with 2 5 < M < 2 20 ) for the critical value fi c = 1 : (a) 
Probability distribution of the maximal weight w max in the region 1/2 < w max < 1 (b) Probability distribution of Y2 in the 
region 1/2 < Y 2 < 1 



6. Conclusion : comparison with the REM and the DPCT 



In this Appendix, we have described the weight statistics in Levy sums for the critical value fi c = 1, to compare 
with the REM and the DPCT cases studied in the text. Although the three models have the same properties in 
the low-temperature phase with /i = T/T c < 1, we find here that the three models have different critical finite-size 
properties. The REM and the Levy sums involve a single finite-size exponent 



vrem = 2 (B15) 

VLevy = 1 (B16) 



and both have a singularity exponent 



&REM — — (JLevy (B17) 

in continuity with its low-temperature value a = fi — 1 — > 0. On the contrary, the DPCT involves two exponents 

vdpct = 2 (B18) 

v'dpct = 1 (B19) 

The exponent v = 2 governs the thermodynamics, in particular the entropy and the specific heat, whereas v' = 1 
governs the Yfe statistics. Moreover, the singularity exponent at criticality 

a DPCT 1.5 (B20) 
is very different from the limit of its low-temperature value a — T/T c — 1 — > 0. 
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